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ABSTRACT 

Using a sample of 92 UV continuum-selected, spectroscopically identified galaxies with (z) = 2.65, 
all of which have been imaged in the Lya line with extremely deep narrow-band imaging, we examine 
galaxy Lya emission profiles to very faint surface brightness limits. The galaxy sample is represen- 
tative of spectroscopic samples of LBGs at similar redshifts in terms of apparent magnitude, UV 
luminosity, inferred extinction, and star formation rate and was assembled without regard to Lya 
emission properties. Approximately 45% (55%) of the galaxy spectra have Lya appearing in net ab- 
sorption (emission), with ~ 20% satisfying commonly used criteria for the identification of "Lyman 
Alpha Emitters" (LAEs) [Wo(Lya) > 20 A]. We use extremely deep stacks of rest-UV continuum and 
continuum-subtracted Lya images to show that all sub-samples exhibit diffuse Lya emission to radii 
of at least 10" (^ 80 physical kpc). The characteristic exponential scale lengths for Lya line emission 
exceed that of the Aq = 1220 A UV continuum light by factors of 5 — 10. The surface brightness 
profiles of Lya emission are strongly suppressed relative to the UV continuum light in the inner few 
kpc, by amounts that are tightly correlated with the galaxies' observed spectral morphology; however, 
all galaxy sub-subsamples, including that of galaxies for which Lya appears in net absorption in the 
spectra, exhibit qualitatively similar diffuse Lya emission halos. Accounting for the extended Lya 
emission halos, which generally would not be detected in the slit spectra of individual objects or with 
typical narrow-band Lya imaging, increases the total Lya flux [and rest equivalent width Wo(Lya)] 
by an average factor of ^ 5, and by a much larger factor for the 80% of LBGs not classified as LAEs. 
We argue that most, if not all, of the observed Lya emission in the diffuse halos originates in the 
galaxy H II regions but is scattered in our direction by H I gas in the galaxy's circum-galactic medium 
(CGM). The overall intensity of Lya halos, but not the surface brightness distribution, is strongly 
correlated with the emission observed in the central ~ 1" - more luminous halos are observed for 
galaxies with stronger central Lya emission. We show that whether or not a galaxy is classified as a 
giant "Lyman a Blob" (LAB) depends sensitively on the Lya surface brightness threshold reached by 
an observation. Accounting for diffuse Lya halos, all LBGs would be LABs if surveys were sensitive 
to 10 times lower Lya surface brightness thresholds; similarly, essentially all LBGs would qualify as 
LAEs. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift 
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Although the Lyman a (Lya) emission line of neu- 
tral H is expected to be prod uced in pr odigio us amounts 
by star-form ing galaxies (e.g. [Partridge k, Peebleslll967l : 
iMeiedfTOTGh , it has long been appreciated that the astro- 
physics affecting observations of Lya are far more com- 
plex than for oth er lines of abundant species due to res- 
onant scattering (lSpitz eilll978t iM eier fc TerlevTchI 119811: 
iCharlot fc Faimi993f) . The very large cross-section in the 
Lya transition means that emission from a gas cloud or 
nebula may have been strongly altered in intensity, kine- 
matics, and apparent spatial distribution by the time 
it reaches an observer. Similarly, information about the 
initial source of observed Lya emission may be lost or ob- 
scured, with the apparent source simply being H I gas re- 
sponsible for scattering in the observer's direction. Con- 
sequently, the dominant process producing Lya emission 
may often be ambiguous; possibilities include photoion- 
ization by young stars or AGN, line emission following 
collisional excitation of H atoms, or simply scattering 
from intervening H I gas that happens to favor the ob- 
server's direction. 

In the absence of dust, the standard expectation for 
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Lya emission produced in H II regio ns for "Case B" (i.e. , 
ioniz ation-bounded) r ecombination (|Brocklehurstl 119711) 
and a 'ChabrieiJ (|2003[ ) stellar initial mass function (IMF) 
for high mass stars is that each solar mass of star for- 
mation produces a Lya luminosity _L(LyQ!) ~ 2.0 x 10'*^ 
ergs s~^. For the same IMF, the far-UV continuum 
light produced per solar mass of SFR near the wave- 
length of Lya has an expected monochromatic lumi- 
nosity in the range 40.0 ^ log L\^cont ~ 40.3 ergs s~^ 
(jLeitherer et al.lll999D . The predicted rest equiv- 
alent width of Lya emission is then given by (Lya) ~ 
L(Ly a)/LA.cont ^ 100 - 200 A (see also ICharlot fc Falll 
119931) , with values near the lower end of this range ex- 
pected for continuous star formation lasting more than 
^ 3 X 10^ yrs, roughly the mini mum dynamical t imescale 
for L* LBGs at - 2 - 3 (e.g.. lErb et alll2006bn . Under 
the above assumptions, the period of time over which 
Lya emission has Woihya) > 100 A would be very 
brief, after which the line-to-continuum ratio reaches 
an asymptotic value of Wo(Lya) — 100 A. Thus, for a 
UV continuum-selected sample, one would expect only 
a small fraction of galaxies to be caught during a time 
when their intrinsic VFo(Lya) exceeds 100 A^^. 

When dust is mixed throughout the scattering 
medium, one expects selective extinction of Lya pho- 
tons compared to those in the nearby UV contin- 
uum due to the much larger effective path length 
traversed by a line photon before escaping in t o the 
intergalactic me dium (e.g.. JMeier fc TerlevichI 119811 : 
iHartmann eFallllQSllNeufeldlllQQOO . This effect is often 
cited when observed Lya emission lines are much weaker 
than the Case B expectations discussed above (e.g.. 
Chariot fc Fa!? 1993: Shap lev et all l200l IHaves et alJ 



strongly affected. 



2010HKornei et al...2010 1. Since most continuum-selected 



high redshift galaxies in current spectroscopic surveys 
appear to have at least some dust, and the vast major- 
ity have Lya equivalent widths Wn (Lya ) < 100 A (e.g., 
iShaplev et all 120031: iKornei et al.ll2010( ). this conclusion 
would seem reasonable. On the other hand, it is also 
possible, at least in principle, for Lya photons to experi- 
ence less attenuation by dust than continuum photons, in 
the case of a clumpy ISM in which dust is located only 
within the clumps which are rarely penetrated by Lya 
photons (Ncufcld 1991; Finkelstein et al. 2008). There is 
no reason to believe that the two competing effects could 
not both be at work within different regions of the same 
galaxy. 

Even without dust, however, resonant scattering pro- 
duces spatial and/or spectral diffusion of Lya photons 
leading to emergent line emission whose properties de- 
pend on the geometry, kinematics, and H I optical 
depth distributions within the gase ous circumgalactic 
medium (CGM) surrounding a galaxy (jSteidel et al.ll2010l 
[S2010]). In the zero-dust case, the total Lya luminos- 
ity would be unaltered by resonant scattering, but, as 
we detail below, the detectability of Lya could be very 

^'^ Note that this value is a factor 1.8 higher than would be 
obtained assuming a Salpeter (1955) IMF because a given number 
of ionizing photons is associated with a smaller total SFR for the 
Chabrier IMF. 

For a sample selected by Lyo (as opposed to continuum) 
emission, this may not be the case. 
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Fig. 1. — Comparison of Lya line equivalent widths measured 
from spectra compared to those inferred from Cont-NB colors in 
deepLya imaging. The imaging measurements use isophotal aper- 
tures defined by the extent of Lya flux to a surface brightness limit 



■ , which is typical of the 



of ~ 1 - 2 X 10" 
deepest Lya narrow-band imaging surveys. 



In S2010, we characterized the distribution of cool 
gas in the CGM of star-forming galaxies with redshifts 
2 ;S z 3 and attempted to understand the kinematics 
and line strength of the ISM absorption and Lya emis- 
sion in the context of galaxy-scale gaseous outflows. In 
brief, we found that UV-selected galaxies within a fac- 
tor of a few of L* in the far-UV continuum luminosity 
function (corresponding at z ^ 2.5 to appa.rent m agni- 
tudes 7^ ~ 24 - 24.5- see iReddv fc Steidell l2009h have 
a CGM that can be traced by H I (Lya and Ly/3 ab- 
sorption) and several strong absorption lines of metallic 
species (e.g., C II, C IV, Si II, Si IV) to galactocentric 
distances of ^ 120 kpc using the spectra of faint back- 
ground galaxies. The measurement used more than 500 
galaxy pairs on angular scales 1 — 15" to map out the 
absorption line strength as a function of galaxy impact 
parameter h (i.e., the physical separation of the two lines 
of sight at the redshift of the foreground galaxy) for each 
observed species. In slit spectra of the CGM "host galax- 
ies" , the bulk of observed Lya emission, when present, 
is almost always strongly redshifted, while the strong in- 
terstellar (IS) absorption lines are strongly blue-shifted. 
S2010 presented a geometric and kinematic model that 
reproduces many of the observed trends. 

In the context of the model, Lya photons escape the 
galaxy in an observer's direction mainly by scattering 
from optically thick H I gas located on the far side of the 
galaxy's stars, but having the same overall (outflowing) 
kinematics as the IS gas seen in blue-shifted absorption. 
We used the transverse information from the galaxy pairs 
combined with line-of-sight information available from 
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Fig. 2. — Comparison of the Lya line equivalent width distri- 
bution from spectroscopic measurements versus that inferred from 
CB-NB colors in Lyo imaging. The imaging measurements use 
colors within isophotal apertures defined by the extent of Lya flux 
to a surface brightness limit of ~ 1 — 2 X 10~^* ergs cm~^ 
arcsec"'^, which is typical of the deepest Lya narrow-band imag- 
ing surveys. The statistics are for the mean and standard deviation 
(left) of individual values (left), and the median and inter-quartile 
range (right) for each set of measurements. 



tfie gafaxies' own far-UV spectra to construct a consis- 
tent geometric and kinematic modei of galaxy scale out- 
flows in the context of a very well-studied population 
of high redshift star-forming galaxies. That is, we com- 
bined the line profiles of IS absorption lines and Lya 
emission in the galaxy spectra themselves (sampling the 
kinematics and line strength for galactocentric impact 
parameter b ^ 0) with IS line strength measurements 
at & >> (using close angular pairs of galaxies) to in- 
fer the 3-dimensional distribution of CGM gas surround- 
ing an average galaxy in the spectroscopic sample. We 
suggested that the CGM gas seen in absorption would 
also constitute a scattering medium through which Lya 
photons must traverse in order to be observed. High 
velocities and large velocity gradients together with gas 
covering fraction /c < 1 through much of the CGM al- 
low Lya photons to diffuse spatially outward, favoring 
escape of Lya photons last scattered (in the observer's 
direction) from atoms with velocities well off resonance 
with respect to any H I that remains between the loca- 
tion of the last scattering and the observer. If true, one 
might then expect to observe scattered Lya emission over 
the same spatial scales for which strong HI and low-ion 
metallic absorption is seen, i.e., ~ 80 — 90 kpc, even if all 
Lya photons originated in the galaxy's H II regions. 

Clearly, scattering will substantially modify both the 
spatial and spectral distribution of Lya photons emer- 
gent in a particular direction, and at the very least may 
cause Lya emitting regions to appear distinct from the 
UV continuum emission even if both share a common 
origin. Slit spectra commonly optimized for the com- 



pact size of the continuum emitting regions of typical 
star-forming high redshift galaxies may encompass only 
a fraction of emergent Lya emission. The relevant angu- 
lar scale for the optically-thick CGM H I gas is ~ 10" 
(~ 80 physical kpc at z ^ 2.5), whereas a typical ex- 
traction aperture for a slit spectrum is ~ 1'.'2 x 1'.'4 - 
a difference of a factor of more than 180 in solid angle. 
Thus, even if the Case B-expected production rate of 
Lya photons were to escape the CGM of a galaxy, it is 
likely that the emission would be distributed over such 
a large region that a narrow slit would miss most of the 
Lya flux; even very deep narrow-band images might leave 
much of the flux unaccounted-for due to limited surface 
brightness sensitivity. 

In this paper, we present direct observational evi- 
dence showing that extended Lya scattering "halos" are 
a generic property of high redshift star-forming galax- 
ies, including those that have no apparent Lya emission 
lines in their far-UV spectra. In §2 we describe a sam- 
ple of 92 UV-continuum-selected galaxies for which both 
rest-far-UV spectra and deep narrow-band Lya images 
are available, and discuss the relationship between Lya 
properties measured using both techniques. In §3 we use 
composite UV spectra, as well as Lya and continuum 
image stacking, to measure Lya emission extending to 
very low surface brightness thresholds for various galaxy 
sub-samples. The results and their implications for the 
nature of Lya emission in star-forming galaxies are de- 
scribed in §4, discussed in §5, and summarized in §6. 

Throughout the paper we assume a Lambda-CDM cos- 
mology with V,m = 0.3, r^A = 0.7, and h = 0.7. 

2. THE GALAXY SAMPLE 

The galaxies used in this paper are drawn from 3 survey 
regions where we used UV-color selection t o select galax- 
ies with 1 . 8 ^ z ^ 3.4 for spec troscopy (jSteidel et al.l 
[20031 [200l iShaplev et all [20051) . In addition to com- 
pleting extensive "Lyman Break Galaxy" (LEG) spec- 
troscopic follow-up, we have also imaged the 3 regions 
using narrow-band (NB) filters centered at the observed 
wavelength of Lya at the redshift of galaxy over-densities 
we had previously identified from the continuum-selected 
spectroscopic sample. Table [1] summarizes the NB obser- 
vations in these fields, all of which are among the deepest 
NB images ever obtained for Lya at redshifts z ~ 2 — 3. 
The number of continuum-selected galaxies with spec- 
troscopic redshifts falling within the redshift range sub- 
tended by the NB filter bandpass in each field are also 
summarized in Table [T] 

We now briefiy comment on each of the fields observed: 

• SSA22a has a redsh if t "spi ke" centered at z = 3.09 
([Steidel et al.l 119981 |2003[ ) which was first iden- 
tified from the spectroscopic follow-up of LBGs. 
It w as first ima g ed in Lya at the same redshift 
by Steidel et al.l (j2000D . who discovered two very 
large (> 100 kpc) "Lya Blobs", prompting sev- 
eral subsequent studies of Lyg-sele cted objects us- 
ing d e eper narrow-band data (e.g.. iMatsuda et ahl 
[25ol IHavashino etall [20p: INestor et al.l 120 111) . 
Here we include the 22 continuum color-s elected 
LBG s with spectroscopic redshifts (Stcidel ~et al.l 
[20031: IShaplev et all [20061 ) lying within a 5'. 5 x 7^6 
region with especially deep Lya NB observations 
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TABLE 1 

Summary of Lyo Narrow-Band Observations 



Field NB^ N<= Telescope/Instrument Date texp"^ PSF<= 5(Lya,obs)' ^(Lya, z = 2.65)^ 

HS1549-I-195 4667/88 2.802 - 2.875 27 Keck 1/LRIS 2010 May 18,000 0'.'86 1.29 1.59 

HS1700+643 4018/90 2.266 - 2.340 43 Palomar 5m/LFC 2007 Jul 80,280 1'.'20 2.30 1.54 

SSA22a'i 4980/80 3.063 -3.129 22 Keck 1/LRIS 2005 Aug 33,880 0'.'80 0.94 1.50 

Subaru/SuprimeCam 2002 Sep 25,800 

^ Central wavelength/bandwidth of NB filter, in A. 
Redshift range included between NB filter half-power points. 

Number of continuum-selected, spectroscopically identified galaxies with NB measurements. 

Total integration time, in seconds. 
° Stellar FWHM in arc seconds after smoothing to match the CB and NB PSF prior to photometry. 
' Observed surface brightness isophotal threshold (1.5(t), in units of 10~^* ergs cm~'^ arcsec"'^. 
s Isophotal surface brightness threshold, corrected to z = 2.65, in units of 10~^* ergs s~^ cm~^ arcsec"'^. 

The field imaged with Keck/LRIS is a 5^5 X 7'.6 subset of the LBG survey field from ISteidel et al.l (12003 . |2000| ). The NB image is a 
combination of the LRIS images and archival Subaru images, discussed by Nestor et al 2011. 



(Table [I]) 

• HS 1700-1-64 is a survey field centered on the the po- 
sition of a hyper- luminous (r' = 16.0, or L ~ 1.1 x 
IQ^'^Ijq) z = 2.751 QSO. A galaxy over-density was 
again identified from spectrosco pic follow-u p, with 
z = 2 .299 ±0.03 (Sha plev et al. 2005; Stcidc leraLl 
12005') ■ We have subsequently obtained very deep 
NB imaging in both Ha and Lya at this redshift 
(Erb et al 2011, in prep.). We include in the 
present sample the 43 continuum-selected galaxies 
with spectroscopic redshifts placing the Lya tran- 
sition within a NB filter designed for follow-up of 
the proto-cluster. 

• HS1549+195 is another survey field centered on 
the position of a hyper-luminous QSO (r' = 15.9, 
or L ~ 1.4 X 10"Lo), with zq = 2.842. Once 
again a galaxy over-density was identified from the 
LBG spectroscopic follow-up, in this case centered 
on the redshift of the QSO itself The NB4670 
filter was designed to follow-up on the galaxy over- 
density, and in response to the serendipitous dis- 
covery of spatially offset (and plausibly fluorescent) 
Lya emission associated with a z — 2.842 damped 
Lya absorption (DLA) absorption system identi- 
fied in the spectrum o f a faint background QSO 
(|Adelberger et al.l [20061 ) . 

In all 3 fields used in the present paper, galaxies were 
selected using rest-UV (LBG) color selection and ob- 
served spectroscopically u sing the Keck 1 1 0m telescope 
and LRIS spectrograph (|Oke et al.l 119951 : ISteidel et al.l 
120041 ) prior to the Lya imaging, so the resulting sam- 
ple should be relatively unbiased with respect to Lya 
properties. The full sample of 92 galaxies with mean 
redshift (z) = 2.65 is broadly representati ve of UV- 
selected spectroscopic samples (^e.g.. .Steidel et al.ll2003l: 
Shaplev et aIll2003HSteidel et al.ll2004HAdelberger et al.l 
2004f ) in terms of both continuum and Lya properties: 
for example, they have 23.4 < TZab ^ 25.5 with median 
(mean) of TZab = 24.47 (24.50), and spectroscopically- 
measured Wo(Lya) in the range —37 A (absorption) 

^® As we will show below, there is no evidence that the presence 
of the QSO has significantly altered the overall Lya emission of the 
galaxies at the same redshift. 



to - 1-89 A (emi ssion) with medi an Wo(Lya) ~ +0.9 A 
(cf. IShaplev et al. 2003; Rcddy e FalllfoOSl : iKornei et al] 
I2OI0I ). The spectroscopic measures of Lya are based on 
extraction apertures of angular size 1'.'2 (the slit width) 
by ~ 1'.'35, independent of wavelength, so that Lya 
and the UV continuum light are measured over identi- 
cal spatial re gions. We used a m ethod similar to that 
described by iKornei et al.l ()2010D to measure Wo(Lya) 
directly from the galaxy spectra. 

Lya equivalent widths (and fluxes) were also measured 
for the same set of 92 galaxies using a comparison of 
deep narrow-band (NB ) and continuum (C B) images. 
As discussed by (e.g.) ISteidel et all (|2000f ). care must 
be taken since spectroscopic and imaging measurements 
of Lya may not be measuring the same quantities. For 
measurements of Lya line emission from CB-NB color, 
the photometric aperture is often defined by the region 
within an isophote corresponding to a particular Lya 
surface brightness threshold, which of course depends on 
the depth of the Lya image. It also depends on the suit- 
ability of the continuum measurement for estimating the 
UV continuum flux density in the vicinity of the Lya line, 
which may require a color-term correction and / or correc- 
tion for Lya line contamination. The 3 Lya images used 
here are comparably deep to the deepest Lya surveys 
to date, with la surface brightness thresholds of 1.53, 
0.86, and 0.63 xlO""'^* ergs s~^ cm~^ arcsec"^ for the 
HS1700, HS1549, and SSA22 fields, respectively^^. Al- 
though these observed SB thresholds differ by a factor of 
more than 2, the deeper data at higher redshift result in 
rest- frame Lya surface brightness thresholds which differ 
by less than 10%. The last column of Table [T] shows the 
relative surface brightness thresholds when all 3 datasets 
are shifted to the mean redshift of (z) — 2.65. 

Important to generating Lya line images for rela- 
tively continuum-bright galaxies (and for measuring line 
equivalent widths independently of spectroscopy) is a 
measure of their far-UV continuum (hereinafter CB, or 
mAB[1220A]) near the wavelength of the Lya line. Ide- 
ally, the CB should have the same effective wavelength 
as the NB without including the Lya line itself. For 
the 3 fields presented here, deep CB images were cre- 
ated using linear combinations of two broad-band fil- 

Isophotal apertures corresponding to 1.5<t above the local sky 
were used for NB-selected catalogs in all 3 cases. 
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ters bracketing tfie NB passband; detaiis of tliese proce- 
dures are described in tlie individual papers cited above. 
Briefly, in HS1700 we used very deep C/„(3550/700) and 
0(4730/1100) images obtained in 2001 May with the 
William Herschel 4. 2m Telescope prime focus imager (see 
iShaplev et al.l 120051 ) to create a "UG" continuum image 
with an effective wavelength of 4010A. For SSA22a, we 
used archival B and V images taken with the 8.2m Sub- 
aru telescope with Suprimecam to create a "BV" CB 
image with Acff = 4980 A (see Nestor et al 2011). The 
HS1549 field was treated somewhat differently, since the 
deepest broad band image (10,800 s integration with 
Keck/LRIS) was obtained in the V band"'^^ using the 
LRIS red channel contemporaneously with the March 
2007 NB4670 images on the blue channel. A less-deep G- 
band image (2500 seconds with Keck/LRIS-B) was used 
to estimate the appropriate (object-dependent) color cor- 
rection needed to adjust the deeper V-band images to 
an effective wavelength near 4670 A. Since the observed 
range in continuum color among the sample galaxies at 
a given redshift is small (e.g., the mean and standard de- 
viation in observed broad-band color for the 43 z ~ 2.30 
galaxies in the HS1700 field is (G - 7^) = 0.26 ± 0.12, 
and {Un — G) = 0.80 ± 0.20 where the standard devia- 
tion in both cases includes photometric scatter), and the 
passbands are separated by only ~ 300 A in the galaxy 
rest frame, we believe that systematic errors associated 
with producing the CB image at the appropriate effective 
wavelength is likely very small (<< 0.1 mag.). 

In brief, the CB and NB images were first scaled to 
have matching zeropoints based on photometry of spec- 
trophotometric standard stars and by calculating the rel- 
ative system throughput in each filter passband as a 
cross-check. The suitably scaled CB or NB images were 
smoothed using a Gaussian kernel to match the stellar 
point-spread functions (PSFs) in the two images; the fi- 
nal PSF size for each field is listed in Table [TJ Matched 
aperture photometr y was performed using dual image 
mode in SExtractor (jBertin fc Arnoutd[l996[) , with CB- 
NB colors measured using photometric apertures defined 
in the NB image at an isophotal threshold equivalent 
to 1.5(7 per pixel above the estimated local sky back- 
ground. The CB zero point was iteratively adjusted by a 
small amount (< 0.1 mag in all cases) so that the median 
color of all objects in the image having 23 < CB < 26 
(the principle range expected for the galaxies of interest) 
has CB-NB=0, corresponding to identical flux density 
measured in each band. The statistical error in the mea- 
surement of CB-NB can be conservatively estimated from 
the dispersion in color for all objects in the same range 
of apparent magnitude, which is quite small because of 
the intrinsically narrow range in color and the depth of 
both the CB and NB images [cr(CB - NB) ~ 0.05 - 0.1 
mag.] A continuum-subtracted Lya line image (here- 
inafter "Lya" image) was formed by subtracting the 
scaled continuum image from the NB image 

Because the data were obtained using a dichroic with a tran- 
sition wavelength of ~ 5000A (dSOO), the V passband was shifted 
to slightly longer wavelength (A^ff ~ 5506 A instead of 5464 A. 

Among the 3 fields, only the SSA22a CB includes a small 
overlap (~ 2.5% of its full bandwidth) with the NB Lya passband; 
this would have the effect of a small (negligible for our purposes) 
over-subtraction of the continuum when producing the Lya line 



We measured VF(Lya) (in units of A) from the CB-NB 
color using the simple relationship 



W^o(Lya) -B 



NB 



-[^g0.4|CB-NB| _ ^ 



CB-NB 
ICB -NBI 
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where Bj^b is the appropriate rest-frame bandwidth in 
A of the NB filter (19.6 A, 23.4 A, and 27.3 A for SSA22, 
HS1549, and HS1700, respectively). Note that with this 
definition, WQ{Lya) — when CB — NB, and positive 
(negative) values indicate net Lya emission (absorption). 
At an isophotal threshold of ~ 1.5 x 10^^* ergs s^^ cm^^ 
arcsec^^, the typical solid angle subtended by the de- 
tection isophote in the narrow-band (NB) image for the 
continuum-selected galaxies is ~ 3.8 arcsec^, ~ 2.3 times 
larger than the spectroscopic aperture. If the spatial dis- 
tribution of Lya emission is significantly different from 
that of the continuum light, then the measured colors 
[and hence the inferred W^o(Lya)] could differ from the 
spectroscopic values. Figure [1] compares the measure- 
ments of Wo(Lya) from the spectra versus those based 
on the CB-NB colors for the same 92 galaxies in the 
current sample, while Figure [2] compares the two distri- 
butions. There is a modest tendency for the value of 
Wo{Lya) measured from the NB images to be larger in 
absolute value (whether in absorption or emission) near 
the extremes of the distribution, though they have very 
similar mean and median values (Figure [2]) and agree 
well when \ Wo\ is small. If significant Lya flux were dis- 
tributed on still larger angular scales (with lower Lya 
surface brightness) while the same is not true of the UV 
continuum light, then even the larger NB-based Wo(Lya) 
would underestimate the true values. 

3. INFERENCES FROM STACKED COMPOSITES 

3.1. Spectroscopic Stacks 

In order to measure Lya emission with SB well be- 
low the detection threshold for individual objects, we 
constructed composite spectra and images after divid- 
ing the sample of 92 into several subsets, summarized 
in Table [2l We used the values of Wq (Lya) measured 
from the CB-NB colors for all galaxies, with apertures 
deflned by the isophotal thresholds listed in column 9 of 
Table [T] This method has generally smaller statistical 
uncertainties and aperture corrections compared to the 
spectroscopic measurements, and facilitates comparison 
with most deep Lya surveys, which are based primar- 
ily on equivalent widths and fluxes inferred from the NB 
photometry 

As illustrated in Figure [21 the median Wo(Lya) 
from both NB imaging and spectroscopic measure- 
ments is close to zero, in agreemen t with previous re- 
sults for continuum-selected samp les (iSteidel et al.|[2000l: 
IShaplev et alll2003l: iKornei et al.ll2010[ ). In forming sub- 
sets of the sample of 92, we used the NB Lya measure- 
ments to split the sample into "Lya Em" , those that 
have Lya in net emission (52), and "Lya Abs", those 
having Lya in net absorption (40)--see the second and 
third rows of Table [2l respectively. Two other subsets 

image. 

'^^ We have verified that none of the results of this paper depends 
significantly on whether the imaging or spectroscopic measures of 
Lya are used to define the subsets. 
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TABLE 2 

LY-a AND Continuum Surface Brightness Profiles for Composites 



Sample^ 


Number 




(1220 A) 


(10-18) 


(kpc) 


(10-18) 


fee" 

(kpc) 


F(Lya)d 
(10-16) 


Ltot{Lya)° 
(10*2 ergs s-l) 


Wo{Lya, specY 
(A) 


Wo{Lya,tot)s 
(A) 


All 


92 


2.65 


24.60 


2.4 


25.2 


87.2 


3.4 


1.7 


9.7 


+ 6.9 


+36.0 


Lya Em 


52 


2.66 


24.40 


3.1 


25.6 


136.3 


2.9 


2.5 


14.3 


+13.2 


+44.9 


Lyo Abs 


40 


2.63 


24.72 


1.5 


20.8 


52.5 


4.5 


0.7 


4.0 


- 4.4 


+16.8 


All non-LAE 


74 


2.65 


24.56 


1.4 


25.5 


124.9 


2.8 


1.4 


8.0 


+ 1.0 


+29.1 


LAE only 


18 


2.64 


24.68 


3.9 


28.4 


110.3 


2.9 


4.0 


22.8 


+29.2 


+92.9 


Lyo Blobs 


11 


2.59 




15.7 


27.6 






11.5 


65.7 







^ Galaxy sub-sample, drawn from the full sample (All) of 92 continuum-selected galaxies with Lya imaging. The details of the sub-samples are described in 
the text. 

^ Average continuum apparent magnitude at Ao — 1220 A, estimated from the CB photometry. 

^ Best fit parameters assuming SB profile S{r) — C'nCxp{ — h/bn), where Cn is in units of IQ-^^ ergs s^^ cm^^ arcsec — ^. The sub-scripts / and c refer to the 

Lya line and UV continuum profiles, respectively. 

^ Average integrated Lya flux, in units of IQ-^^ ergs s" ^ cm-^. 

^ Average integrated Lya luminosity, in units of 10^^ ergs s"!, assuming {z) — 2.65. 
1 Lya rest equivalent width measured from spectrum (Figure [31. 



^ Lya rest equivalent width of total Lya flux, 

were made consisting of galaxies satisfying the criteria 
commonly adopted for "Lyman a Emitters" (LAEs), i.e., 
Wo (Lya) > 20 A, of which there are 18 (20% of the 
total), with the remainder (74 of 92, or 80%) placed 
in a sub-sample called "non-LAEs", i.e., all continuum- 
selected LBGs that would not be selected as LAEs. 

For each sub-sample listed in Tabled a composite far- 
UV spectrum was created by shifting the observed, flux- 
calibrated spectra into the galaxy rest-frame using the 
prescriptions given in S2010, re-sampling the rest-frame 
spectra to 0.5 A pix"^ , and averaging. Each stacked com- 
posite spectrum was scaled so that the continuum level 
near Lya matched that obtained from the photometric 
stack of the same subset of galaxies, discussed below. 
The correction was typically a factor of ~ 1.5, and was 
applied for the sole purpose of placing the continuum 
levels for the same subsets on the same flux scale. The 
resulting stacked spectra are shown in Figure [31 mea- 
sured properties of the composite spectra are given in 
the figure and listed in Table [2] 

3.2. Lya and CB Stacks 

A 25" x 25" sub-image ("postage stamp") centered 
on the position of the continuum centroid of each galaxy 
was extracted from the CB image and the (continuum- 
subtracted) Lya images after scaling them to a com- 
mon zero point as discussed above. Masks were created 
by performing object de tection on each continu um sub- 
image using SExtractor (jBertin fc Arnoutsll 19961 ). These 
were used to exclude pixels lying within the detection 
isophotes of any object other than the central one, and 
were applied during the stacking to both the CB and Lya 
images. Two stacked images were formed for each sub- 
set listed in Table [5] (straight averages, with masking), 
one for the CB image and another for the Lya line im- 
age. Figure m compares the average CB image with the 
Lya image for the full sample of 92 galaxies, while Fig- 
ure [5] compares the azimuthally-averaged surface bright- 
ness profiles of the same composite CB and Lya images. 
Figure m shows clearly that, on average, Lya emission is 
detected to radii of at least 10" , or ~ 80 physical kpc 
at (z) = 2.65. Figure [5] shows that the average CB light 
profile for the same galaxies is much more compact and 
drops below the SB detection threshold for 6 ^ 20 kpc 



( ^ 2"5). Figure [5] also shows what the Lya line profile 
would look like if Wo(Lya) = lOOA (i.e., Case B) and 
Lya and CB light had the same spatial distribution on 
average. 

The surface brightness profiles for both the continuum 
and the Lya line images are reasonably well-fit by an 
exponential of the form S{b) — Ciexp{~b/bi) for pro- 
jected radii beyond the central arcsec; the parameters 
of the best-fit values for the normalization Ci and scale 
length hi are given in Table [2] for each sub-sample as 
well. In the full image stack, the effective surface bright- 
ness detection thresholds are a factor of ~ 10 lower than 
for individual galaxies; it is clear that the distribution of 
Lya emission is very different from that of the continuum 
for every sub-sample, with best-fit Lya scale lengths of 
bi ~ 20 — 30 kpc compared to the corresponding con- 
tinuum emission which has 6c — 3 — 4 kpc. It is im- 
portant to note that the true difference in scale length 
is larger, since we have made no attempt to de-convolve 
the profiles from the seeing disk, which was FWIIM~ 
0.86, 1.20, and 0.80" for HS1549, HS1700, and SSA22a, 
respectively. The continuum profiles of the stacked com- 
posite CB images have FWHM ~ 1"2 - 1"4, indicating 
average (seeing de-convolved) galaxy continuum sizes of 
FWHM~ 0.80" (ct ~ 0.35" ). These de-convolved angu- 
lar sizes are also consistent with measureme nts of simi- 
lar galaxies in deep H ST/ACS images (e.g.. iPeter et al.l 
[20071: iLawlFall [20071 1 

The stacked Lya and CB images as in Figures [5] and 
[Hi represent unweighted averages of all galaxies in the 
sample (with masking as described above). This choice 
was motivated by the desire to preserve the photometric 
integrity of the stacks so that fluxes could be measured 
directly using aperture photometry, but also because any 
scaling or weighting would require deciding what the rel- 
evant figure of merit should be. Medians are often used 
to suppress outliers in stacked data sets, but they have 
the disadvantage for the present application of not pre- 
serving flux in two-dimensional images, of working best 
when scaling has been applied to individual images go- 
ing into the stack, and of suppressing real signal as it 
approaches the noise level. Nevertheless, in Figure [3 we 
show a comparison of the surface brightness profiles for 
median-combined stacks as compared to mean-combined 
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Fig. 3. — Composite spectra formed from the average within the subsamples detailed in Table [2] Within each panel, the number of 
galaxies going into the stack is listed after the sub-sample name; the second line in the annotation lists the rest-frame equivalent v^idth of 
the Lya line measured from the composite spectrum, with the convention that positive values indicate net emission. 



for both line and continuum. 

We have argued that our sample of galaxies with 
(z) = 2.65 has emission line and continuum properties 
characteristic of those in the full LBG spectroscopic sur- 
veys at these redshifts. Figure [5] shows that the diffuse 
Lya emission is also consistent among the 3 survey fields 
taken individually. This is important, since each field 
samples galaxies at a different redshift with observations 
subject to a different set of conditions, using different 



telescope and instrument combinations. We also note 
that the bright QSO known to lie within the survey vol- 
ume in the HS1549 field appears not to have had a sig- 
nificant effect on the Lya emission from the galaxies in 
our sample^ ^. 

The galaxy regions most likely to be affected by excess ionizing 
radiation from the QSO would lie in the outer parts; whether this 
radiation would increase or decrease the amount of hya emission 
from galaxies would depend on the physical state of the gas. The 
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Fig. 4. — (Loft:) Scaled far-UV continuum image produced (as described in the text) from tfie average of 92 continuum-selected LBGs, 
drawn from 3 independent fields. The regions shown are 20" (~ 160 kpc physical at 2 = 2.65) on a side, with a grid spacing of 2" . 
(Right:) The continuum-subtracted, stacked hya image for the same sample of galaxies. In both panels, the contours are logarithmically 
spaced in surface brightness with the lowest contour shown at ~ 2.5 X 10~^^ ergs cm~^ arcsec"^. 



Figureinishows the measured surface brightness profiles 
and best-fit exponential parameters for the sub-sample 
with net Lya emission (see the second row of Table [2]). 
The profile of the LyaEm composite is qualitatively sim- 
ilar to that of the full sample. The main difference is in 
the central Lya surface brightness — Lya Em objects 
have (on average) Lya surface brightness well above the 
threshold for individual detection (light shaded region in 
Figure. [5] and [6]) to projected distances of 6 ^ 20 — 25 
kpc (2'.'5— S'.'l), whereas even the peak Lya SB in the 
full sample (Figure [5]) barely reaches the (individual) de- 
tectability threshold. Also illustrated in Figure [6] is the 
SB profile for the "Lya Abs" sub-sample. The Lya SB 
scale length for the Abs sample is still ~ 4 times larger 
than for the corresponding continuum light (third row 
of Table [2])- and has a comparable Lya scale length bi 
to that measured for the Lya Em sample (second row 
of Table [2]) , albeit with a ~ 3 times lower normalization 
for & ^ 20 kpc. Clearly the difference is much larger for 
6 1" , where the Abs sample exhibits a large "hole" in 
which Lya absorption strongly dominates. 

Figure [9] reproduces the SB profiles of the "Lya Em" 
and "Lya Abs" sub-samples together with the average 
profile of LAEs (green) , and non-LAEs (cyan) . The LAE 
sub-sample is very similar to that of the larger "Lya Em" 
subset, but has an average SB a factor of ~ f .5 higher for 
6 XL 10 kpc (and a factor ~ 3 higher for 6 < 10 kpc.) Also 
plotted for comparison (Figure El) is t he average SB pro- 
file of 11 "Lya Blobs " (LABs; see e.g. ISteidel et al. l l2000l : 
iMatsuda et al.ll20M ). which for the present we define as 
Lya selected objects with detected isophotal diameters 
d > 5" , discovered in the same 3 survey fields. None 
of the 11 Blobs is included in the main galaxy sample, 
since they do not have central continuum sources that 
satisfy the usual LBG color criteria. As indicated in Ta- 
ble El the Blobs have an average Lya luminosity ~ 7 

galaxies in the HS1549 field are fainter by about 30% on average 
(in terms of apparent continuum magnitude) than in the other two 
fields; the differences in the Lya profiles on small scales may be a 
consequence of this selection issue. 
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Fig. 5. — The observed average surface brightness profile for 
the 1220 A continuum light (blue) and the Lya line (red) for the 
full sample of 92 continuum-selected galaxies, evaluated over the 
same rest-frame bandwidth sampled by the Lya image (24.3 A). 
Note that these profiles are simply the azimuthal averages of the 
stacked images shown in Figure [4| The light-shaded region indi- 
cates the range of typical Lya surface brightness threshold reached 
by deep Lya surveys for the detection of individual objects. The 
dashed lines show the surface brightness profile assuming that 
S{b) = C;exp(— fe/6;) with parameters given in Table [2] The cor- 
responding angular scale at (z) = 2.65 is given along the top axis. 
For the purpose of comparison, we also show the Lyo profile ex- 
pected for the same sources under the assumption of "Case B" Lya 
to CB ratio, no destruction of Lya by dust, and no spatial diffusion 
of Lya photons due to resonant scattering (i.e., the Lya and CB 
profiles would be identical in shape, and since H'o(Lya) = 100 A, 
the Lya line image would be a factor of ~ 4.1 brighter than the 
continuum in the effective rest-frame bandwidth of 24.3 A.) 
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Fig. 6. — Same as FigureO comparing the average surface bright- 
ness profiles for the the sample divided according to whether the 
NB measurements indicate net Lya "Abs" or "Em" . 
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Fig. 7. — Same as FigureO comparing the surface brightness pro- 
file measured from median-combined rather than mean-combined 
stacks of Lya and CB sub-images. Note that the line-to-continuum 
ratio is not necessarily preserved in the median stacks. The gen- 
eral effect of this alternative processing is to decrease the measured 
Lyo scale lengths by ~ 25%, though the Lyo profiles remain much 
more extended than the continuum profiles. 

times higlier tlian an average galaxy in our sample. It 
appears tliat even tlie most extreme LABs do not have 
fundamentally different SB profiles compared to those 
of typical galaxies in the sample except that their sur- 
face brightness normalization exceeds the typical detec- 



tion threshold to 6 ~ 50 kpc ('^ 6" ). In other words, 
if one were routinely sensitive to a surface brightness of 
^ IQ-w Qj-gg g-1 cm~^ arcsec~^ , all continuum- selected 
LBGs would he "Blobs"'^'^. We will return to a more de- 
tailed discussion of Lya demographics in §4. 

If we assume for the moment that the extended Lya ha- 
los represent photons originating in the galaxies' H II re- 
gions, we can use the composite Lya line and CB images 
to measure the integrated Lya line-to-continuum ratio, 
usually parametrized as Wo(Lya), the Lya equivalent 
width. The total Lya fluxes have been measured directly 
from the calibrated stacked images (column 8 of Tabled]); 
a comparison with the continuum flux density measured 
near the wavelength of Lya from the CB images (col- 
umn 3 of Tabled allows the calculation of Wo(Lya, tot) 
(column 12). These numbers can be compared directly 
with the spectroscopic measurements (column 11) for the 
same galaxy sub-samples. The values that include the 
diffuse Lya extending to ^ 80 kpc radii around galaxies 
exceed the spectroscopically-inferred II^o (Lya, spec) by 
an average factor of ~ 5 for the full galaxy sample. Fig- 
ure[Tn]shows the cumulative fraction of the total Lya flux 
as a function of aperture radius b (in arc seconds) for the 
galaxy samples in Table [H It is interesting to note that 
including the spatially extended Lya emission brings the 
average galaxy into the range that would nominally qual- 
ify as a LAE {Wo(Lya) > 20 A)- even for the "No LAE" 
sub-sample that explicitly excludes the 18 conventional 
LAEs (row 4 of Table O. 
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Fig. 8. — Average observed Lya surface brightness profiles for 
galaxy sub-samples separated by field, as indicated. The Lya pro- 
files were scaled according to the relative continuum flux density 
in each field, for display purposes. The small differences in the 
mapping of angular scale to physical scale have also been removed 
to facilitate the comparison. The measured values of the average 
Lya rest equivalent width are (Wo(Lya)) = 28.1 A 28.6 A, and 
42.0 A for HS1549, HS1700, and SSA22, respectively. 

■^^ Conversely, if the LABs were several times less luminous 
but had the same surface brightness profile, they would fail to be 
recognized as "blobs" at all. 
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One can also compare the measured large-aperture 
Wo{l-iya) with expectations for the Lya to continuum 
ratio for Case-B recombination and no dust. As dis- 
cussed above, an asymptotic value of Wo{Lya) ~ 100 A 
is expected when star formation has been continuous for 
tsF ~ 4 X 10^ years^^. Note that approximately the same 
value of Wo (Lya) is expected as long as Lya photons do 
not suffer greater attenuation by dust than continuum 
photons just off the Lya resonance. Thus, the fact that 
most of the values of fesc,rei = Wo(Lya)/100A (TableEJ 
are significantly smaller than unity means that Lya pho- 
tons suffer greater extinction than the continuum, by fac- 
tors ranging from 1.1 — 6.0 with an average ~ 2.8 for the 
full sample of 92 galaxies. 

We note that the values in Table El for fesc,rei have 
been obtained using a method that appears to differ from 
that used in some r ecent work (e.g.. [ Gronwall et al.|[2007l : 
iNilsson et all 120091: iKornei et al .''2010'). Most estimates 
of fesc,rei usc stcUar population synthesis models to esti- 
mate the level of extinction, which is then used to derive 
SFR to calculate the expected Lya luminosity based on 
the assumption of Case B recombination and the form 
of the stellar IMF. While we are using largely identi- 
cal SED modeling to estimate SFRijv,corr, we use the 
observed Wo(Lya) as a direct observational estimate of 
fesc,rei- The present method relies on the same assump- 
tions about the stellar IMF (i.e., based on a Salpeter-like 
IMF for high mass stars) to estimate the Lya photon 
production rate per unit star formation; the difference is 
that we rely on the ratio of Lya photon production to 
that of Aq — 1220 A continuum photons from the same 
ensemble of stars. The advantage of using the equiva- 
lent width measurement is that it should be independent 
of extinction if Lya and ~ 1220 A continuum photons 
experience the same attenuation, and would directly re- 
flect the relative attenuation of line and continuum if 
^Lya 7^ ^1220 ■ Bccause both of these methods rely on 
measuring the integrated Lya line flux, an underestimate 
of Lya relative to the continuum will cause fesc,rei to be 
underestimated by the same factor. For the present pur- 
poses, we prefer using the method relying on 11/o(Lya) 
since it depends on a largely independent measurement 
that may avoid propagating possibly large systematic er- 
rors in the estimates of E(B-V) from SED fitting [or from 
the assumed extinction curve, which relates E(B-V) to 
A(A)] to the calculation of fesc,rei- In general, we expect 
that the under-counting of Lya photons due to the aper- 
ture effects discussed above are likely to dominate any 
differences in inferred fesc,rei- 

Column 10 of Table [2] compiles the average Lya lu- 
minosity Ltot(Lya) implied by the measured value of 
-Ptot(Lya) assuming the sample mean redshift (z) = 2.65. 
If one naively converts these numbers to an equivalent 
star formation rate (i.e., divide by 2 x 10**^ ergs s~^ to 
yield SFR in units of Mq yr~^) the results range from 
2.1 — 11.8 M0 yr~^ for the various sub-samples, with an 
average of 3.1 M© yr~^ for the full sample. In Table[3]we 
have compiled the statistics of the far-UV inferred SFRs 

The typical inferred age fo r galaxies similar to those in the 
prese nt sample is ~ 500 Myr fe.g.. lShaDlev et al.l2005l : IReddv et al.l 



and continuum extinction for each of the sub-samples 
from Table [21 The extinction estimates are parame trized 
by E{B - V) and assume the ICalzetti et al.l (|2000) star- 
burst attenuation curve; E(B-V) was estimated from 
SED fits when available, and using the far-UV contin- 
uum slope for the ~20% of galaxies lacking adequate 
near-IR photometric coverage for SED fitting. In the 
context of the assumed starburst attenuation relation, 
the extinction (in magnitudes ) at 1500A is A(1500A) — 
ll.lQE(B- V) (Meurcr erZIIT99l[Caketti et al.l[2?)0fll : 
IReddv et al.) I2b06. ,201011 The UV continuum mag- 
nitudes near rest-fra me 1500 A were used to esti- 
mate SFRnv (e.g., iM adau et al.' '1998'; iSteidel et all 
119991 lAdelberger fc Steid el 2000) with a median value 
of {SFRuv) — 6.3 Mq yr~^. Applying the median 
E(B-V) to the median SFRjjv within each sub-sample 
implies that 11.3 < {SFRuv,corr) ~ 46.5 Mq yr-\ with 
an overall median of SFRijv,corr = 34.3 Mq yr~^ — 
very cl ose to the mean o f the LEG sample observed in 
Ha by lErb et al.l (|2006b[ ) and consistent with the mean 
bolometric luminosity of identically-sele cted LBGs esti- 
mate d using multiple SFR indicators (Re ddv et al.ll2006l 
|2010[) . The median E{B — V) varies considerably among 
the sub-samples, so that the median attenuation of the 
UV continuum is inferred to range from ~ 2.5 for the 
LAEs to ~ 7 for the Lya Abs sub-sample. In Table [3] 
we list the inverse of this factor, which we have called 
fesciUV) = SFRuv/SFRuv,coir- Our estimate of the 
fraction of all Lya photons produced by photoionization 
in the galaxy H II regions that have been detected is 
then given by /esc,tot(Lya) = fosc,roi(Lya) x fcsc(UV). 
These values range from /esc, tot (Lya) ~ 0.37 for the 
LAE sub-sample to /esc, tot (Lya) ~ 0.024 for the "Lya 
Abs" sub-sample. The average for the entire sample 
is /esc,tot(Lya) ~ 0.061. We note that this fraction 
is close to the ave rage value of /esc (Lya) estimated by 
iHaves et al.l (|2010H based on a very different approach in- 
volving a comparison of Ha and Lya luminosity density 
at z ~ 2.2. 

The last column of Table |3| shows the inferred ra- 
tio j4(Lya)/E(B — V), where both quantities are ex- 
pressed in magnitudes and E(B — V) is inferred from 
the stellar SED. The numeric value of this ratio is ~ 
19.5 ± 1.5 for all sub-samples except the LAEs, which 
have A(Lya)/E(B - V) ~ 11.9. Since the ICalzetti et al.l 
(I2OOO,) extinction curve predicts that the UV continuum 
near the Lya line has A{1220A) / E{B -V) ~ 12, it seems 
that Lya emission from the LAEs drawn from our LBG 
sample exhibit no evidence for selective extinction of Lya 
photons, while for other LBGs the attenuation A{Lya) 
is ~ 1.6 times higher than for continuum photons for 
the same value of E{B — V). If the Lya escape frac- 
tion is controlled by processes confined to H II regions, 
the result suggests that E{B — V)neb — V E{B — V)stars 
with T] ~ 1.6, on average. This can be compared with 
the relationship inferred for nearby star-forming galaxies, 
r/ ~ 2.5, based on measurements of the Balmer decre- 
ment (jCalzetti et al. 2000). At present, there are few 
galaxies for which Lya Ha, and H/3 have all been mea- 
sured, though there are some indications that the same 
value of E{B — V) applies for both continuum starlight 
and Ha for galaxie s similar to those in the current sample 
(|Erb et al.ir2006at but see iForstcr Schreiber et al, .20091 
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Fig. 9. — As for Figure [S] where here the Lya radial SB profiles 
are shown for each of the sub-samples in Table [2] along with the 
exponential models for each. Also included for comparison is the 
average surface brightness profile of 11 giant Lya "Blobs" observed 
in the same 3 survey fields (red). 

for possibly conflicting evidence). 

4. IMPLICATIONS OF DIFFUSE Lya HALOS 

Diffuse Lya emission from the outer parts of actively 
star forming galaxies is an unavoidable consequence of a 
gaseous CGM so long as some component of it is optically 
thick to Lya photons and some fraction of Lya photons 
initially produced in H II regions are not absorbed by 
dust at smaller galactocentric radii. Calculation of the 
emergent Lya emission is undoubtedly complex, since 
it will depend on the details of the gas-phase structure 
and kinematics as well as the relative distribution of the 
sources (e.g., H II regions) and the sinks (e.g., dust) of 
Lya photons. A fully successful model requires 3-D ra- 
diative transfer calculations and all of the relevant spatial 
and kinematic information as input. Such a treatment is 
far beyond the scope of this paper; however, it is inter- 
esting to ask whether the spatial profiles of Lya emission 
from the same star-forming galaxies can be understood 
in the context of a schematic model. In this section, we 
describe such a model that begins with inferences on the 
structure and kinematics of CGM gas from S2010, and 
then test for consistency with both the Lya emission ob- 
servations and the absorption-based S2010 CGM model. 

4.1. A Model for Lya Scattering Halos 

We first consider the probability that a Lya photon 
produced in a galaxy's central few kpc will escape in the 
direction of a particular observer's line of sight. The es- 
cape probability will depend on the kinematics and opti- 
cal depth distribution of the CGM gas, and so one might 
expect it to be closely related to the characteristics of 
absorption lines observable both in the galaxy spectra 
themselves {b ;$ 2-3 kpc) and in lines of sight to back- 
ground objects with at larger impact parameter 6 >> 0. 




2 4 6 8 10 



b (arcsec) 

Fig. 10. — The cumulative fraction of the large-aperture Lyo 
fiux as a function of angular aperture radius b for each galaxy sub- 
sample. The vertical dashed line drawn at 6 = 0'.'65 indicates the 
typical effective aperture for the slit spectra of the same objects. 
The dark blue curve corresponds to the cumulative continuum fiux 
(for the stack of the full sample) as a comparison. 

The conditions necessary for a Lya line photon to es- 
cape in the direction of a particular observer are: 1) it 
must either be emitted at a frequency that is well off res- 
onance for any H I in the foreground (i.e., between the 
point of emission and the observer), and/or 2) it must be 
scattered in a direction that happens to have low spatial 
covering fraction fc of H I^''. 

For extended Lya produced by scattering in a gaseous 
halo, the observed surface brightness profile S{b) will 
then be related to the integral along the line of sight 
at impact parameter b of the product of a) the Lya pho- 
ton density, b) the probability that a Lya photon will 
be scattered in our (the observer's) direction, and c) the 
probability that once scattered in our direction a pho- 
ton will proceed to escape the nebula before being scat- 
tered once again. The situation is somewhat analogous 
to the galaxy outflow model used to match absorption 
line equivalent widths Wq vs. impact parameter b pre- 
sented in S2010. Figure [TT] shows the assumed geometry 
(cf. Figure 23 of S2010.) In the absorption line 
of sight to a background object pierces the radial flow 
at projected distance 5, and the resulting absorption line 
strength is modulated by the integral along the line of 
sight of the quantity 1 — /c(r, Vout), where r is the galac- 
tocentric radius and Vout (r) is the flow velocity at radius 
r. As discussed by S2010, the velocity field in the ab- 
sorbing gas can have a large effect on the strength of ab- 
sorption lines in the spectra of background sources when 
the transition is saturated, even if the covering fraction is 
significantly smaller than unity. S2010 argued that con- 

^■^ In the limit of no H I gas outside of a galaxy's H II regions, 
the emergent Lya line would have roughly the same spatial extent 
as that of the UV continuum. 
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TABLE 3 

Inferred Continuum and Lyq Attenuation for Sub-Samples 



Sample^ 




E(B - Y)" 


SFRuv'^ 


SFRuv.corr*^ 


foBc(UV)f 


fcsc,rol(Lya)^ 


fcsc,tot(Lya)'' 


A(Lyo)/E(B - V)' 


All 


24.47/24.50 ± 0.55 


0.17 


6.3 


34.3 


0.17 


0.36 


0.061 


17.9 


hya Em 


24.50/24.55 ± 0.56 


0.11 


6.0 


18.6 


0.32 


0.45 


0.144 


19.1 


hya Abs 


24.44/24.43 ± 0.54 


0.19 


6.6 


46.5 


0.14 


0.17 


0.024 


21.3 


All non-LAE 


24.43/24.42 ± 0.54 


0.18 


6.7 


42.6 


0.16 


0.29 


0.046 


18.5 


LAE only 


24.85/24.80 ± 0.50 


0.09 


4.5 


11.3 


0.40 


0.93 


0.372 


11.9 



^ Galaxy sub-sample, drawn from the full sample (All) of 92 continuum-selected galaxies with Lya imaging. The details of the sub-samples arc 
described in the text. 

^ Median and mean/standard deviation of continuum apparent magnitude at Aq — 1500 A 
Median E(B - V) inferred from SED fitting. 

d 



Median SFR, in Mq yr , from UV continuum with no dust correction. 
^ Median SFR after correction based on E(B-V) and Calzetti (2000) reddening relation. 
^ Fraction of 1500 A photons escaping galaxy. 

^ Relative escape fraction of Lya photons, = Wq (Lya)tot /lOO A. 
^ Fraction of hya photons escaping, fcsc.tot — fcsc,rcl X fcsc(UV). 

^ Ratio of attenuation of Lya photons to E(B — V) when both are expressed in magnitudes. 

sistency between the absorption line strength as a func- 
tion of impact parameter on one hand, and the strength 
and profile shape of lines observed in the spectra of the 
galaxies themselves on the other, requires large veloci- 
ties and velocity gradients in the gas. The absorption 
cross-section is dominated by outflowing material, and 
the flows are inferred to be clumpy (i.e., multi-phase), 
with both high- and low-ionization ionic species observed 
over similar ranges of velocity and galactocentric dis- 
tance. In the context of the CGM model, most of the 
acceleration of cool gas to high velocity occurs in the 
inner several kpc. The covering fraction fc of gas giv- 
ing rise to absorption in a particular transition decreases 
with increasing galactocentric distance r, modeled as a 
power law of the form fc oc r~'^ . 

In general, the larger the range of gas-phase bulk ve- 
locity sampled along the observer's line of sight at im- 
pact parameter b, the greater the chance that a scattered 
Lya photon will reach the observer without further scat- 
tering. For simplicity, in our model we assume that all 
scattering events are isotropic, and that the gas-phase ve- 
locity field is axisymmetric and is a monotonic function 
of galactocentric distance r (see S2010 for a justification 
of this assumption). If the bulk velocity field in the out- 
flow has a range and amplitude much larger than that 
of the local velocity dispersion in the H I gas^^, or when 
when fc << 1, the problem can be reduced to a geo- 
metric one in which the covering fraction depends only 
on galactocentric radius r, fc = fc{f)- Clearly it would 
be interesting to measure the velocity field of extended 
Lya emission in order to gauge the role kinematics play 
in the transfer of Lya photons. Unfortunately, beyond 
the central, high surface brightness regions there are few 
constraints on the line shapes, and at present we have 
only (projected) spatial information integrated over the 
full range of velocity. 

In the models discussed here, this is assumed to be the case 
based on the results presented in S2010. While velocity is not 
used explicitly as a model parameter, large velocity gradients along 
the line of sight directly affect the probability that a scattered 
Lya photon will ultimately escape. In other words, the effective 
covering fraction of optically thick H I as seen by a Lyo photon 
emitted from a particular position in the CGM implicitly includes 
an integral ove r velocity even if it is not explicitly used as a model 
parameter. See lSteidel et al1l2010l for a more detailed discussion of 
this issue. 




LOS 



Fig. 11. — Coordinate system for the schematic model of scat- 
tered Lya emission from galaxies, where b is the line of sight im- 
pact parameter, r is the galactocentric distance, and ^?e// is the 
assumed physical size of the scattering medium. In the model, 
Lya photons are produced near r ~ 0, after which they diffuse 
outward until they are either destroyed or they escape the scatter- 
ing medium. The covering fraction of gas at galactocentric radius 
r is assumed to be of the form fdr) oc r~"i (see also S2010). 



When considering emission (rather than absorption in 
the spectra of background objects) one needs to account 
for the Lya "source function" which varies with spatial 
position, as well as variations in opacity parametrized by 
/c(r). The Lya photon density available to contribute 
to the observed S'LyQ(^) will depend on the fraction of 
Lya photons that have been able to diffuse outward to 
r ^ b, which may be only a small fraction of the Lya 
photons initially produced by recombination in H II re- 
gions. When the covering fraction is high at small radii, 
one would expect the emergent Lya emission from that 
region to be suppressed - photons are either destroyed 
or radiatively trapped until they make their way to lo- 
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cations from which escape is more probable. The flux 
of Lyof photons (assumed to be produced at small r at 
a ^ constant rate related to the SFR) at galactocentric 
radius r will be reduced by an overall geometric factor 
1 / 47rr^ , and by the destruction of Lya via absorption by 
dust grains^^ . 

This diminished Lya radiation field would produce no 
Lya halo at & > r if /c(r) = 0, since photons would 
appear to be released from a Lya "photosphere" that 
would approximately extend to the edge of the gas dis- 
tribution. The apparent outer edge of the Lya scattering 
halo should correspond to the radius at which /c becomes 
negligible and the scattered component of Lya falls be- 
low the observational threshold. At small radii, where 
/c(r) — 1 and the optical depth encountered in any direc- 
tion is substantial, (e.g., for galaxies having absorption- 
dominated continuum spectra), Lya photons will be res- 
onantly trapped for a large number of scattering events 
before diffusing spatially outward. Most of the dust ab- 
sorption, if present, would be expected to occur in such 
regions. Once fc{T) falls below unity at larger radii, Lya 
photons which have not been destroyed may be scattered 
in the observer's direction. Thus, the relative rate of Lya 
scattering events at radius r will be cx fc{r)/ATTr'^ where 
/c(?') is the H I covering fraction. The chance that a 
scattered Lya photon will be emitted in the observer's 
direction (without any further interactions prior to es- 
cape) increases with decreasing characteristic /c, with 
probability roughly cx [1 — fdr)] for fc{r) < 1. The Lya 
surface brightness as seen by an observer in a particular 
direction will then be proportional to the product of these 
two terms, integrated along the line of sight through the 
galaxy at impact parameter b: 

Q ,M o /■+'— /,(r)[l-/,(r)] 

JLya[b) (X So J dl (2) 

where I is the coordinate distance along the observer's 
line of sight at impact parameter b, Imax — {R^s — b'^)'^ , 
5*0 is a normalization for the surface brightness distribu- 
tion, and i?eff is the effective size of the scattering halo 
(see Figure [TT]) . Note that the integrand tends toward 
zero when fc{r) ~ 1, qualitatively accounting for the 
suppression of Lya emission in regions with high fcif). 
Clearly, when /c = 1 the purely geometric model is no 
longer valid, since the Lya emission intensity associated 
with an optically thick region can be "negative", i.e., 
there is a net removal of Lya photons at that spatial po- 
sition that will reduce the net surface brightness along 
that particular line of sight. 

In the S2010 CGM model, the radial dependence of the 
covering fraction of gas was found to be consistent with a 
power law of the form fc{r) oc , where 0.3 ^ 7 0.6 
depending on the ionization state of the tracer ion; fc{r) 
becomes consistent with zero for r > Res ~ 90 kpc for 
most of the observed ions. If we take the power law form 
for fc (with index 7) characteristic of the highest H I 
optical depth material, equation [2] can be used to predict 
Shyaib) given an overall normalization S'o, a character- 
istic radius tq at which fc{r) first falls below unity [i.e. 

Under the assumption of spherical symmetry, Lya photons 
scattered at smaller radii are returned to the "pool" of Lya photons 
potentially available for scattering at larger radii. 



fc{i") = (r/ro)"'''], and the the effective size of the scat- 
tering medium, Reg. We account qualitatively for the 
variation of the spatial profile of Lya emission in the cen- 
tral regions of a galaxy by (artificially) allowing /c > 1 
at r < To, using an extrapolation of the same power law 
form for /c(f). Since the integrand becomes negative 
when /c > 1, this leads to suppression of the Lya sur- 
face brightness for any line of sight that intercepts such 
a region. In practice, the central Lya emission must be 
substantially suppressed to match the observed profiles 
of any of the galaxy subsets, including the LAEs (cf. 
Figure O . 

Figures[T2]and[T3]show example models based on equa- 
tion [5] compared with the observed composite Lya sur- 
face brightness profiles; the corresponding model param- 
eters are given in Table ID The overall shape of the pre- 
dicted SB profile is sensitive to the value of 7 parametriz- 
ing the radial dependence of the covering fraction of H I; 
7 < 0.3 produces Lya profiles that are fiatter than ob- 
served, while 7 > 0.8 predicts Lya emission which falls 
too rapidly with increasing b fFigure [T5|) . As discussed 
above, the shape of the central portion of Si^yaib) is mod- 
ulated by adjusting the galactocentric radius ro where 
fcifo) = 1 (i-e., To serves as a normalization of the maxi- 
mum covering fraction) . The presence of a central "hole" 
in Si^ya{b) (as observed for the Lya Abs sub-sample in 
Figure I12p can be reproduced by increasing tq so that 
the transition from net Lya absorption to net Lya emis- 
sion moves to larger galactocentric radius. Once Lya 
has a finite probability of escape (i.e., where fdr) < 1 in 
the context of our simple model), the residual Lya pho- 
tons at r > ro become available for re-direction toward 
an external observer who then perceives the photon be- 
ing "emitted" from a position at impact parameter b in 
projection. 

The values of 7 required to produce model profiles in 
reasonable agreement with their observed counterparts 
(Tablc[4]and Figure [T2|) are near the high end of the range 
inferred from the behavior of absorption line strength Wq 
versus impact parameter b (S2010). One possible expla- 
nation for slightly steeper profiles is that the emission 
models assume no Lya photons are destroyed once they 
propagate beyond r ~ ro; if Lya has a finite chance of 
being absorbed by dust at r > ro, the additional attenu- 
ation of Lya would manifest itself as a steepening of the 
profile with respect to the pure scattering model. That 
the Lya Abs model exhibits both the steepest decline in 
Lya surface brightness (7 ~ 0.8 compared to 7 ~ 0.6 for 
the other sub-samples) and the largest global extinction 
correction (§3 and Tabic [3]) suggests dust may not be 
confined solely to the central regions in such galaxies. 

The model of the CGM proposed by S2010 almost cer- 
tainly does not provide a unique explanation for the IS 
absorption line strength and kinematics as observed in 
the spectra of background galaxies; however, we have 
shown, with a simple extension of the model, that scat- 
tering of Lya photons from the same CGM gas can can 
also account for Lya emission with radial surface bright- 
ness profiles and physical extent consistent with the ob- 
servations. Regardless of the model details (which admit- 
tedly could be incorrect), the very similar physical scales 
involved (i?cff — 90 kpc) suggest a close causal connec- 
tion between the cool gas observed to produce strong H I 
and low-ionization metallic absorption lines in the spec- 
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Fig. 12. — Same as Figure O but where models described by 
equation [2] have replaced the exponential profile used in Figure |9] 
Parameters for the 4 models shown are summarized in Table |4] 
All 4 model curves have Reft = 90 kpc and fc{r) oc (r/ro) ' 
with 7 ~ 0.6 — 0.8. The "Lyo Abs" model produces a central 
hole in the Lyo emission by adjusting the normalization of fc with 
the parameter ro; the larger value of ro indicates that the CGM 
remains optically thick to Lya photons to larger galactocentric 
radii than for the other sub-samples. 



TABLE 4 

Parameters for Model Lya Spatial Profiles 



Sample^' So"" ro (kpc)<= 7^ R^g (kpc)' 



All 


11.5 


2.2 


0.6 


90 


Lya Em 


17.0 


2.0 


0.6 


90 


Ly« Abs 


4.5 


5.9 


0.8 


90 


All non-LAE 


7.0 


2.0 


0.6 


90 


LAE only 


25.0 


1.9 


0.6 


90 



^ Galaxy sub-sample, drawn from the full sample (All) 
of 92 continuum-selected galaxies with Lyct imaging. 
^ Intensity normalization for model (see Eq. [2j, in units 
of 10""^* ergs s~^ cm"^ Hz~^ 
^ Galactocentric radius at which /c — 1 

Power law index in the radial behavior of the covering 
fraction, /c — (r /rQ)~"' . 

^ Effective size of CGM region producing detectable Lya 
emission, in kpc. 

tra of background continuum sources, and spatially ex- 
tended Lya emission from the same host galaxies. 

4.2. Comparison with Lya Emission in Simulations 

A scenario in which extended Lya emission around 
galaxies is dominated by scattering of Lya photons ini- 
tially produced inside the galaxies, rather than by ex- 
ternal processes, has been the focus of a number of re- 
cent galaxy mode ls including treatment of Lya radia- 
tive transfer fe.g.. Verhamme et al."2008': 'Laurse n et al] 
Il009a,b; Zheng et al. 2010a^b; Barnes et al. 2011). Each 
of these studies places emphasis on different aspects of 
the model galaxies, and cursory examination suggests a 
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Fig. 13. — Same as Figure [S] where the preferred model is drawn 
with the solid black curve, corresponding to So = 11.5 X 10~^* 
ergs cm~^ arcsec"^, R^g = 90 kpc, ro = 2.2 kpc, and 7 = 0.6, 
for covering fraction parametrized as fc{r) oc (r/ro) ' (see also 
S2010). The dashed curve shows a model with the same parame- 
ters, except R^ff = 50 kpc, while the dotted curve assumes 7 = 0.2, 
with all other parameters as for the preferred model. 

qualitative similarity to the observations presented here, 
since Lya scattering leads to the spatial re-distribution 
of the Lya emission as seen by an observer. Gas-phase 
kinematics play a large role in determining how much 
Lya emission esca pes the galaxies, and in all of the mod- 
els except those of I Verhamme et all 120081 (which do not 
explicitly consider the spatial distribution of Lya emis- 
sion) the dominant velocity field is associated with in- 
fall/accretion. The predicted Lya line profiles tend to 
be asymmetric and sometimes double-peaked, usually 
dominated by photons that are blue-shifted with respect 
to the galaxy systemic velocity- a configu ration that is 
very rarely observed in galaxy spectra (e.g..'Pe ttini et aP 
I2000t Shaplev et al. 2003; Steidel et al. 201^^ Also, 
while the 3-D models all produce Lya emission that is 
significantly more extended than the UV continuum, the 
predicted surface brightness profiles of scattered emis- 
sion generally declines much more rapidly than for the 
observed LBG Lya halos (i.e., most woul d fall well be- 
low th e current surface brightness limit). B arnes et al.l 
(|201lD have pointed out that higher outfiow velocities 
tend to produce more extended Lya emission in the con- 
text of their models, which include both infiows and out- 
flows of gas, so that perhaps the missing ingredient is the 
presence of higher-velocity outflows than have generally 
been modeled. 

Most of the simulations work on Lya emission from 
galaxies has not highlighted the potential utility of us- 
ing Lya emission observations as a means of revealing 
gas-phase structure in the surrounding CGM and IGM. 
One exception is a a series of recent papers exploring 
how very sensitive NB observations of Lya emission can 
be used along with cosmological simulations (including 
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detailed radiative transfer) to tra ce the under lying large 
sc ale structure at high redshifts ([Zheng et al . 2010a. b). 
In iZheng et al.l (|2010bD . the authors explicitly calculate 
the expected properties of diffuse Lya emission around 
star-forming galaxies, with principal focus on LAEs at 
z ~ 6. Like the scenario we have described above, the 
models assume that the ultimate source of the Lya pho- 
tons seen in emission is the galaxy HIT regions, with ex- 
tended emission resul ting from the d etails of the Lya 
radiative transfer. Zheng et alj ()2010 bl predict that the 
surface brightness profile surrounding individual galax- 
ies will have two distinct components related closely to 
1) the "halo exclusion scale" within comoving distances 
of 0.3h^^ Mpc {b ~ 60 physical kpc at z ~ 6), and 
a larger-scale component arising from galaxy clustering, 
extending to ~ 3 Mpc (comoving), or ~ 400 physical kpc 
at z ~ 6. The smaller scale is similar to the virial radius 
of the characteristic dark matter halos being considered 
in the simulation. It is not completely straightfor- 
ward to move the predictions to z ~ 2.65 for comparison 
with our observations, but (as discussed in S2010) the 
CGM scattering medium observed around z ^ 2 — 3 
galaxies (which we have argued is responsible for ab- 
sorption against background sources as well as for the 
extent of scattered Lya emission) also has a size similar 
to the virial radius r^, ~ 80 — 100 kpc. However, it is 
not clear that the observations are consistent with the 
predictions when it comes to the dependence of the Lya 
emission hal o on oth er properties of the galaxies. In the 
IZheng et al.1 (|2010bD models, the primary driver of the 
surface brightness profile is the gas-phase kinematics of 
the CGM gas; the characteristic scale of the inner com- 
ponent of Lya emission relates to the "infall" region for 
the halo, within which gas is accreting onto the central 
galaxy (the simulations do not have outflowing material, 
and it is the kinematics of infalling material that modu- 
late the escape of Lya photons). On the other hand, in 
our picture the characteristic scale is related to the radial 
dependence of the covering fraction of neutral material 
and the gas-phase kinematics (assumed to be dominated 
by outfiows). Within our sample of LBGs, the stacked 
Lya images of various subsets indicate a rather consis- 
tent exponential scale length of bi ~ 25 kpc with at most 
a weak dependence on Lya or UV luminosity or on the 
fraction of Lya photons that escape the galaxy ISM. If 
there is a trend, it is in the direction opposite to that 
expected in the models. 

Interestingly, many recent theoretical investigations 
focusing primarily on diffuse and extended Lya from 
the outer parts of galaxies or LAEs have deliberately 
neglected the scattering of Lya from the inside out 
(IDiikstra et al.ll2006HDiikstra fc Lo"ebll2009HKeres et al.l 
I2009t IFaucher-Giguere et al.l 120101 IGoerdt et al.l 120091 ). 
Instead, attention has been drawn to Lya emission asso- 
ciated with gas cooling as it accretes onto galaxies ( "Lya 
Cooling" ) , or on Lya fiuorescence as a means of measur- 
ing the intensity of sources of ionizing photons at high 
redshifts. Both of these processes are discussed in §5 

Note that our procedure of masking out all identified con- 
tinuum sources other than the central one when producing the 
L^ct and continuum stacks (§3) would suppress what Zheng et al. 
1201 Oh ) call the "two-halo term" due to clustering, so our observed 
Lya profiles should be compared only with the "one-halo" , central 
component. 



below. 

In any case, there is no doubt that radiative transfer 
calculations will be key to a full understanding of diffuse 
Lya emission from galaxies. However, it is essential that 
the CGM gas distribution and kinematics in the simula- 
tions match real galaxies. Without the correct gas-phase 
model, even the most sophisticated treatment of radia- 
tive transfer cannot yield a realistic result. The observa- 
tions suggest that possibly important ingredients include 
a CGM that is clumpy on small scales and which has very 
large (non-gravitational) velocity gradients dominated by 
galaxy-scale outfiows. 

5. DISCUSSION 

We have shown above that, on average, LBGs with far- 
UV luminosities 0.3 < (L/L'^jy) < 3 at (z) = 2.65 exhibit 
spatially extended Lya emission to physical radii of at 
least 80 kpc (10" ), even when Lya appears only in ab- 
sorption for regions coincident with the UV continuum 
starlight. Figures [Sj [H and [T2] show that the profiles of 
the Lya emission are quite similar in shape independent 
of the spectral morphology, with the main difference be- 
ing the overall intensity normalization and the presence 
or absence of emission spatially coincident with the con- 
tinuum light (i.e., the inner ±5 kpc). The observations 
suggest that the Lya-scattering CGM may be statisti- 
cally universal, with the main variable being the fraction 
of Lya photons able to emerge from the inner few kpc 
region without being destroyed. For example, the differ- 
ence between the Lya Em and Lya Abs (see Table ^ 
spectrally classified subsets is an overall factor of ^ 5 
in the Lya surface brightness at the full continuum ex- 
tent (Figure ini), beyond which the ratio of 5(6) for the 
two sub-samples remains essentially constant. The scale 
lengths for Lya emission {bi ~ 25 ± 3 kpc) are consis- 
tent among the statistically distinct galaxy sub-samples 
in spite of the fact that the integrated line-to-continuum 
ratio varies by large factors among the same sub-samples. 

5.1. Previous Results on Statistical Lya Detections 

Lya emission with physical extent larger than that 
of a galaxy's continuum starlight is not a sur- 
prising result from a theor e tical perspective ( e.g., 
iBarnes fc Haehneltl [2001 [2OIOI : iLaursen et all [2009a . b) , 
and has been observed and noted in many individ- 
ual cases both in the nearby (e.g. Mas-Hesse et al.l 
[200I iHaves et al.l [20071: [Qsthn et al.l [2001) and high 



redshift (e.g iFranx et'aLl Il997l: iMoller fc W arren lll99S. 
jSteidel et all 120001: iFynbo et al.l 120031: iMatsuda et al l 
120041 : 1 AdelbCT ger et al .1 120061: IOuchi'ei~al1 l2008') universe. 
However, relatively few surveys at high redshift have 
reached adequate Lya surface brightness limits to al- 
low the detection of the very low surface brightness 
levels discussed above. An exception is the extremely 
deep spectroscopic survey for Lya emission conducted by 
Rauch ct al. (2008) [R08]. Using a Lya -selected sample 
distributed over the redshift range 2.7 < z < 3.8, these 
authors noted that extended Lya emission was a com- 
mon feature of the LAEs discovered in their survey. A 
spatial stack of all of the Lya emitting sources exhibited 
significant emission (with threshold ~ 1.5 x 10""'^^ ergs 
s^^ cm~^ arcsec"^) to an angular scale of ^ 4" , or ~ 30 
kpc projected physical radius. The R08 sample, as the 
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authors themselves point out, covers a different range 
of UV luminosity compared to most continuum-selected 
LEG spectroscopic surveys- only one of 27 objects has 

V < 25.5, while 80% our sample (which has a median 

V ~ 25.0) has V < 25.5, although there is is a ten- 
dency for the faintest objects to be among those with 
the strongest Lya emission lines (see Tables [2] and [3]). 

Nevertheless, the average surface brightness profile 
for the R08 Lya-selected sample is remarkably similar 
to that of our continuum-selected sample (e.g., compare 
Figure [5] to Figure 20 of ROB). For objects in our "Lya 
Em" sub-sample (Table [2]), the peak Lya SBs are some- 
what higher than for the R08 sample, while the angular 
extent (at the same limiting SB of 1 x 10""'^^ ergs s~^ 
cm~^ arcsec"^) is ~ 2.5 — 3 times larger in the present 
LEG sample. Within our sample there is a significant 
dependence of Wq (Lya) on apparent UV continuum lu- 
minosity, but the average Lya profiles arc similar, as 
shown in Figure 1141 
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Fig. 14. — A comparison of the continuum and Lya surface 
brightness profiles of the full sample divided into two at the median 
continuum apparent magnitude. The "UV bright" sample is a fac- 
tor of ~ 2.0 times brighter in the continuum than that of the "UV 
faint" sample (CB(Bright)= 24.22 versus CB(Faint)= 24.95), but 
the average hya flux for the UV bright sub-sample is 10% smaller 
than that of the UV-faint sub-sample, i.e. Wo(Lya, bright) = 22.0 
A, while Wo(LyQ, faint) = 48.5 A). 

In any case, it is worth pointing out that, under the 
hypothesis that Lya scattering, and not fluorescence, is 
the dominant process producing the observed Lya halos, 
the scattering medium need not be optically thick in the 
H I Lyman continuum. This means that it is not nec- 
essarily correct to associate the observed physical extent 
of Lya emission with regions having N(HI) > 3 x 10^* 

Moving our continuum-selected sample to the somewhat 
higher median redshift of R08 would result in ~ 50% of our sample 
having V > 25.5. 



cm^^ as R08 have suggested - in principle, N(HI) could 
be 1000 times lower and still remain optically thick to 
Lya photons. 

Perhaps more directly analogous to the results of 
the present sample is the narrow-band Lya survey of 
IHavashino et al.l (j2004D . These authors used deep NB 
Lya images in the SSA22 field, and stacked the Lya im- 
ages of 2 2 z = 3.09 continuum-selected LEGs from the 
survev of lSteidel et al. (2003), of which 19 are m common 
with our current SSA22 sample^^. Indeed, Hayashino et 
al showed that significant emission extends to angular 
scales of at least 4" and that the "ring" in the range 
2-4" often contains as much or more Lya flux than the 
inner 9 < 2 " region. They also stated (but did not 
show) that a stack of the 13 galaxies which did not indi- 
vidually exhibit extended Lya emission results in a sig- 
nificant detection on the same 2-4" scales. Although the 
authors did not discuss what physical mechanism might 
have been responsible for their observation, these results 
clearly provided an early indication of the nature of Lya 
emission in L* galaxies, borne out by our larger and more 
sensitive sample. 

5.2. Has the Whole Iceberg Been Detected? 

The level of sensitivity to low-SE Lya emission at high 
redshifts is unlikely to improve by large factors using 
the current generation of ground-based telescopes, and 
so a natural question would be: How much more is there 
at still lower SB? Many Lya surveys (e.g.. iRauch et al.l 
120081: iBunker et aLlll998[ ) have been designed to detect 
Lya fluorescence induced by the metagalactic radiation 
field at redshifts 2 ^ z ^ 3. The radiation field intensity 
is usually expressed as J^/ ~ 2 — 10 x 10~^^ ergs s~^ 
cm~^ Hz~^ sr~^, where the quoted range indicates the 
dispersion amo ng published observa t ional or theoretical 
estimates ( e.g., 'Shaplev et al.' '2 0061: iBolton et al.l [20051: 
IScott et al. 2000; Fauchcr-Giguc re et al.ll2008[) . The ex- 
pected maximum fluorescent signal at z ~ 2.5 — 3.0 is in 
the range 0.2 — 2 x 10^^^ ergs s~^ cm~^ arcsec"^ if the 
only source of io nizing photons is the general UV back- 
ground (see e.g. iCantalupo et al.l 120051 : iKollmeier et al.l 
120101: iFaucher-Giguere et al.l |2010DT These expecta- 
tions clearly lie at or below the current SB thresh- 
olds of any survey completed to date. The dif- 
ficulty of detecting the fluorescent signal from the 
metagalactic UV fleld has instead inspired several 
searches for fluorescence ne ar bright sources of ioniz- 
ing photons, such a s QS Os (iFrancis & Bland-Hawthorn] 
2004; Cantalupo c faLl [200 5: Adclbcrgcr et al. 200C 
iHennawi et al.ll2009() . The results from such studies have 
been mixed. 

A different argument can be used to suggest that 
fluorescence from the UV background will always be 
overwhelmed by Lya scattering from the CGM of star- 
forming galaxies, at least at 2; 2 — 3. This assertion fol- 
lows from the fact that S2010 found that the total absorp- 
tion cross-section contributed by the CGM of LBGs (us- 
ing Roff = 80 — 90 kpc for the detection of low-ionization 

The new NB image used in the present sample includes both 
archival Subaru data as well as an additional 10 hours' integration 
using LRIS on the Keck 1 telescope, and so is substantially deeper 
(~ factor of 2-3 ) , but covers a much smaller area, than that of 
[Havashino etaU 1 I2OOI) . 
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absorption species) can account for a large fraction of all 
gas with N(H I) ^ 2x 10^^ cm~^ (i.e., r > 1 in the Lyman 
continuum, also known as "Lyman Limit Systems"). In 
other words, any gas of sufficiently high N{HI) to pro- 
duce a detectable signal from fluorescence also lies within 
~ 90 kpc of a star-forming galaxy with properties sim- 
ilar to those in our sample. We have shown that these 
galaxies generically exhibit diffuse Lya emission on the 
same physical scales when a surface brightness threshold 
of S'(Lya) ~ 1 x 10~^^ ergs s~^ cm~^ arcsec"^ is reached. 
Unless the fluorescent Lya signal lies at the very top of 
the allowed range, it will have much lower SB than the 
signal we have attributed to scattering from the inside of 
the galaxy out. 

It is more difficult to assess what fraction of ob- 
served Lya emission may be due to cooling processes 
such as t hose described by a n umber of recent au- 
thors (e.g., 'Diikstra fc Loebl 2009t iKollmeier et al.l [2010l: 
iFaucher-Gi gucrc ct al. 201d: iGoerdt et al.l 120091 ) In 
particular, the predictions of the emergent Lya emis- 
sion from cooling gas accreting onto galaxies are ex- 
tremely sensitive to gas tempera ture ()Kollmeier et al.l 
120101: iFaucher-Giguere et al.ll2010() and to the small-scale 
structure in the gas. As a result, the range in Lya flux 
and SB, as well as the galaxy mass dependence and spa- 
tial distribution of cooling emission, must be regarded 
as uncertain by a factor of ^ 10, with an upper bound 
(based on energetic arguments) that can be as large as 
L(Lya) ^ 10'''' ergs s~', but which under different as- 
sumptions could be as small as ~ 5 x lO'" ergs s~' for a 
galax y with Mhaio — 9 x 10" Mq ([Faucher-Giguere et al.l 
I2010D . approximately the mean halo mass of the galax- 
ies in the present sample (see lAdelberger et al.l 120051 : 
IConrov et al1[200l ISteidel et al.ll2010DT 

The observations appear to argue against a signifi- 
cant contribution of cooling radiation to the detected 
Lya halos, at least on average. We have shown that 
the shape of the observed radial surface brightness dis- 
tribution among the LBGs in the sample is remarkably 
consistent beyond the inner ~ 10 kpc, within which the 
Lya intensity for a given continuum luminosity varies 
by orders of magnitude. Moreover, the overall intensity 
scaling for the Lya emission at large radii is strongly cor- 
related with the behavior of Lya emission in the inner 
5-10 kpc region — at the same continuum luminosity, 
Lya absorption-dominated galaxies (on average) exhibit 
diffuse Lya emission with a factor of 3-4 lower normal- 
ization than galaxies with spectroscopically detected Lya 
emission. In the context of Lya cooling radiation, one 
might expect the extended Lya emission to be strongly 
correlated with gala xy mass and/or SFR since it is be- 
lieved by some (e.g., IGoerdt et al.ll2009() that the bary- 
onic accretion rate ultimately controls the SFR. In this 
scenario, the central region of Lya emission might be 
suppressed by higher H I column densities mixed with 
dust, but the outer regions would have no obvious way 
to "know about" the number of Lya photons being pro- 
duced at smaller radii. Instead, one might expect that 
the brightest Lya halos would be associated with the 
"Lya Abs" sub-sample, since these have a median SFR 
nearly 3 (4.5) times larger than the "Lya Em" (LAE) 
sub-samples. Clearly, the observations are inconsistent 
with this expectation. If on the other hand most or all of 



the Lya emission at all radii originates in the central re- 
gions and is subsequently scattered by the CGM gas, the 
density of photons available for scattering at (for exam- 
ple) r = 50 kpc will be very tightly linked to the number 
of Lya photons that successfully diffuse past r "o kpc, 
beyond which the Lya halos appear "self-similar" . The 
emergent Lya luminosities are entirely consistent with 
the observed level of star formation in the galaxies, and 
are more attenuated than the UV continuum, for all sub- 
samples except the LAEs. It is not necessary to invoke 
sources of Lya emission other than scattering (from the 
inside outward) to account for both the Lya luminosity 
and its spatial distribution. 

Under the scattering hypothesis, and further assuming 
that the scattering medium is self-similar for all galaxies, 
then sub-samples with more luminous Lya halos should 
provide information on the degree to which even the cur- 
rent SB threshold might lead to an underestimate of the 
total Lya flux emergent from a galaxy. To increase the 
dynamic range for detecting diffuse Lya emission, one 
might use the observed properties of giant LABs (which 
are well-detected in the stack to 6 ~ 15" ) to estimate 
how much additional Lya flux may lie beyond the SB 
detection threshold near b ^ 8" for more typical galax- 
ies. Under the assumption that diffuse emission from 
LABs and LBGs has a similar origin and differs only 
in total Lya luminosity, the curve-of-growth for LABs 
(Figure [T0|) suggests that an aperture of radius c± 8 " 
would underestimate the total Lya flux by only ~ 10%. 
Thus, further aperture corrections to the integrated Lya 
would probably leave the values of Wo (Lya) (Table [2]) 
and /esc, re; (Lya) (Table [3]) more or less unchanged. At 
least at z ~ 2.65, the current SB limit appears to be 
sufficient to detect most of the "iceberg" . 

Finally, we note that the differences in the intensity 
of the large-scale diffuse emission among sub-samples 
divided according to their spectral morphology suggest 
that galaxy viewing angle is relatively unimportant (on 
average) for Lya emission; that is, most galaxies are not 
LAEs in some directions but strong Lya Abs systems in 
others, consistent with the inference of generally axisym- 
metric CGM gas distributions inferred from the absorp- 
tion line studies (S2010). 

5.3. IS Absorption, Lya Emission, and the CGM 

Perhaps the strongest corre lation (first explored in 
detail by iShaplev et al.l 120031 for galaxies at z ~ 3) 
among the observed spectral properties of LBGs is be- 
tween the strength of low-ionization IS absorption lines 
and the spectral morphology and equivalent width of 
Lya. Galaxies with the strongest Lya emission (among 
the continuum-selected samples) invariably have much 
wea ker than averag e low-ionization IS absorption lines 
fsee lErb et al.|[20Tol for a well-observed example), while 
those with the Lya appearing strongly in absorption 
have correspondingly strong IS absorption features, often 
reaching zer o intensity over som e or most of the line pro- 
file (see e.g. iPettini et al.ll2003 ) indicating unity covering 
fraction. These trends are easy to understand in the con- 
text of the CGM model discussed by S2010 and extended 
in this paper to cover the expectations for Lya scattering 
and its effects on the observability of Lya emission: both 
the IS absorption lines and Lya line strengths and mor- 
phologies are controlled by the kinematics and geometry 
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of the same interstellar and circum-galactic gas. 

Dust certainly plays a role in determining the fraction 
of both Lya and continuum photons that will end up 
reaching an observer. However, the gas-phase geometry 
and kinematics are more directly responsible for the ob- 
served line strength (and line-to-continuum ratios) in the 
spectra. If a galaxy has strong Lya emission emerging 
from the same region as the UV continuum, it must have 
shallow IS absorption lines; if it did not, then at least 
the spatial distribution of Lya (if not also its integrated 
flux) would be substantially modified it would become 
more spatially diffuse. When a slit spectrum (generally 
a small-aperture measurement) shows very strong and 
deep low-ionization IS absorption lines, including Lya, 
it must be the case that any Lya seen in emission will 
have escaped either from a region spatially distinct from 
the continuum (the subject of this paper), or by way of 
scattering from very high velocity material (see S2010). 
Lya emission seen in spectra which also show strong IS 
absorption will be primarily in the latter category, hence 
the nearly universal systemic redshift of Lya emission in 
LBG spectra. We have emphasized above that any Lya 
photons that are not destroyed by dust will eventually 
find their way out of their host galaxy- but will be much 
harder to detect by the time they do. 

The point is that IS absorption and Lya emission are 
causally intertwined through their mutual dependence on 
the structure and kinematics of the CGM on scales from 
a few kpc to ~ 100 kpc. 

6. SUMMARY 

We have presented observations of a sample of 92 
continuum-selected LBGs at {z) = 2.65 having both 
rest-UV spectra and very deep narrow-band Lya images. 
The sample, which is representative of ~ L* LBGs at 

2 ^ z ^ 3, was used to examine the nature of diffuse Lya 
emission from star-forming galaxies as function of their 
spectral morphologies and NB-inf erred Lya fluxes. By 
stacking both UV continuum and Lya line images for 
subsets of the galaxy sample, we are able to study the 
spatial distribution of Lya and continuum emission to 
much lower surface brightness thresholds (~ 1 x 10""'^^ 
ergs cm~^ arcsec"^) than would be possible for in- 
dividual galaxies. We find: 

1. Relatively luminous star- forming galaxies gencri- 
cally exhibit low-surface brightness Lya emission to pro- 
jected radii of at least 80 physical kpc (^ 10 arcsec). 
The extended emission is present even for the stacks of 
LBGs that would be classified based on their spectra as 
having Lya in net absorption. The Lya line to UV con- 
tinuum ratio is always strongly suppressed in the central 
regions of galaxies relative to "Case B" expectations, but 
beyond galactoccntric radii of r ~ 5 kpc, where the con- 
tinuum light falls off very rapidly, Lya emission begins 
to dominate. 

2. The Lya emitting regions have characteristic expo- 
nential scale lengths 5-10 times larger than the the corre- 
sponding UV continuum emission from the same galax- 
ies. It appears that on average all classes of star-forming 
galaxies in the observed range of luminosity would be 
classified as "Lyman a Blobs" if the observations were 
sufficiently sensitive. Similarly, nearly all galaxies would 
also be classified as "Lya Emitters" if their total Lya flux 
were measured using a sufficiently large aperture. Spec- 



troscopic measurements (or, typically deep Lya narrow- 
band surveys) of Lya emission underestimate the to- 
tal Lya flux [or, equivalently, the rest equivalent width 
Wo(Lya)] by an average factor of 5, and a factor of > 3 
even for those classified as LAEs. 

3. The surface brightness distribution, total flux, and 
scale lengths for Lya emission are all consistent with a 
picture in which most or all detectable Lya emission 
is produced in H II regions spatially coincident with 
the galaxies' UV continuum emission. The Lya surface 
brightness is then modified by scattering from the surface 
of H I clouds that are being driven to large galactoccn- 
tric radii by galaxy-scale outflows. The spatial extent 
of observed Lya emission is then dictated by the spatial 
extent of circum-galactic gas with sufficiently large cov- 
ering fraction to have a finite chance of scattering Lya 
photons in the direction of an observer. 

4. The inferred attenuation of Lya emission from 
continuum-selected LBGs is consistently larger than that 
of the UV continuum (A(Lya) ~ 1.6A(UV)) for all sub- 
samples except LAEs, which have A{Lya.) ^ A(UV). 
Most of the attenuation of Lya emission appears to oc- 
cur within ^ 5 kpc of the continuum centroid of a galaxy. 
While the fraction of a galaxies' total Lya photon pro- 
duction that is able to diffuse beyond ~ 5 kpc varies sub- 
stantially, the scattering halo of cool material at larger 
radii leads to self-similar diffuse Lya halos. A simple 
scattering model for Lya emission was presented, based 
on the structure of the CGM gas inferred from mea- 
surements of absorption in lines of sight passing within 
b < 125 kpc of an ensemble of similar galaxies. The 
model successfully accounts for both the typical size and 
the surface brightness profile of Lya emission. 

5. We argue that scattering of Lya photons from 
circum-galactic gas can account for all of the observations 
of continuum-selected star-forming galaxies, and that the 
observed correlations of the intensity of diffuse Lya ha- 
los with the spectral morphology of the central galaxies 
argues against a significant contribution from Lya cool- 
ing of accreting gas. We also argue that Lya scattering 
processes will always dominate over fluorescence (caused 
either by the mctagalactic ionizing radiation field or by 
ionizing photons from inside the galaxy) in producing 
spatially extended Lya emission. 

6. Lya emission and interstellar absorption line 
strengths arc causally intertwined through their mutual 
dependence on the structure and kinematics of CGM gas. 
Galaxies with strong and centrally-peaked Lya emission 
are expected to be associated with shallow IS absorption 
lines, while strong Lya absorption lines that completely 
absorb the UV continuum light of the host galaxy force 
Lya emission to larger galactoccntric radii before escap- 
ing the host galaxy. 

The detection of diffuse Lya emission halos at the 
current surface brightness level has required the equiv- 
alent of ~ 200 — 1000 hours' integration time with a 
lOm-class telescope (accounting for the effective integra- 
tion time of the stacked Lya line images). Given the 
(1 -|- z)^^ dependence of observed Lya surface bright- 
ness (for a given physical luminosity surface density), it 
is not feasible at present to obtain similar results at sig- 
nificantly higher redshifts. However, vastly increasing 
the number of continuum-selected LBGs with sensitive 
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NB Lya observations could in principle trace hya emis- 
sion from galaxy halos until they become indistinguish- 
able from the background. Together with observations of 
the cool gas phase via absorption lines in the spectra of 
backgromid sources (both galaxies and QSOs), such Lya 
emission observations considerably enhance our ability 
to observe directly the distribution of cool baryons and 
their flow rate into and out of forming galaxies during an 
undoubtedly crucial (but not well-understood) period in 
cosmic history. 
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